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ABSTRACT 

Early time-series photometry from NASA's Kepler spacecraft has revealed a planet transiting 
the star we term Kepler-4, at RA = 19'^02'"27".68, S = +50°08'08".7. The planet has an orbital 
period of 3.213 days and shows transits with a relative depth of 0.87 x 10~'^ and a duration of 
about 3.95 hours. Radial velocity measurements from the Keck HIRES spectrograph show a reflex 
Doppler signal of 9.3li[ g ms~^, consistent with a low-eccentricity orbit with the phase expected 
from the transits. Various tests show no evidence for any companion star near enough to affect 
the light curve or the radial velocities for this system. From a transit-based estimate of the host 
star's mean density, combined with analysis of high-resolution spectra, we infer that the host star 



is near turnoff from the main sequence, with estimated mass and radius of 1.223 



1.487 



+0.071 
-0.084 



+0.053 
-0.091 



Mgand 



Rq. We estimate the planet mass and radius to be {Mp,i?p} = { 24.5 ± 3.8 Mq, 
3.99 ± 0.21 i?©}. The planet's density is near 1.9 gcm~^; it is thus slightly denser and more 
massive than Neptune, but about the same size. 
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1. Introduction 

Transiting extrasolar planets provide unpar- 
alleled opportunities for detailed study of the 
physical characteristics of distant solar systems. 
Since the first transiting planet detection a decade 
ago, ground-base d surveys such as TrES, HAT, 



and S uper- WASP (lAlonso et al.ll2004l : iBakos et al 



20041 : IPollacco et all l2006h and the s paceborne 
telescope CoRoT (jBaglin et al.l 120071 ) have lo- 
cated more than 50 transiting planets, spanning a 
large range in size and mass. With the advent of 
NASA's Kepler Mission, we have a new and ex- 
traordinarily sensitive tool for studying transiting 
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planets. Kepler has enough sensitivity to detect 
Earth-size planets orbiting in the habitable zones 
of Su n-like stars during i ts planned 3.5-ycar mis- 
sion (|Koch et al.|[2O10l : iBorucki et al.ll2010l ). The 
first science data to return after Kepler's launch 
in Mar 2009 were time series from a 9.7-day com- 
missioning run, followed after a 1.6-day gap by a 
33.5-day science run. Here we describe the tran- 
siting planet Kepler-4b, one of several transiting 
planets discovered during these first two observing 
intervals. 

2. Observations, Analysis, Tests for False 
Positives 

Observations of the Kepler target field com- 
menced 1 May 2009; the data that we describe 
here are long cadence (LC) photometry, which 
correspond to integration times of 29.426 minutes. 
For a full description of the Kepler field of view, 
obs erving modes, and d a ta processing p i peline . 



obs ervmg modes, and a a ta processmg p i peline , 
see 'Jenk ins et all (|2010l ): ICaldwell et ail (|2010l ). 
For Kepler target stars brighter than r = 13, the 
RMS photometric precision attained for relative 
flux time series is typically better than 2 x 10~ ^ 
per 29-minute integr ation (|Gilliland et al.ll2010l l. 



We detrended photometry from the mission data 
reduction pipeline and searched it for significant 
trai isit-like events using the proced ures described 
by Jenkins et al. ( 2010l ) and by iBatalha et al 
()20m 



One of the transiting planet candidates iden- 
tified by the process just described was the star 
we now term Kepler-4. This star is uncommonly 
bright by Kepler standards. Its Kepler magnitude 
(AB magnitude averaged over the Kepler band- 
pass) is Kepmag = 12.211. Characteristics of this 
star are given in Table 1; briefly, it is a somewhat 
metal-rich {[Fe/H] = +0.17 ±0.06) star of nearly 
solar temperature (5857 ± 120 K), evidently seen 
near the end of its main-sequence lifetime, as ex- 
plained below. 

Figure 1 shows the light curve for Kepler-4, 
folded with a period of 3.21346 dQ Since the 
transit signal was clear, we proceed ed with follow- 
up observations as described in iGautier et al 



^Time series of the photometry and of ra- 
dial velocity data presented here may be re- 
trieved from the MAST/HLSP data archive at 

[littp://archive.stsci.edu/prepds/kepler_hlsp| 
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Fig. 1. — The phased light curve of Kepler-4b 
containing 13 transits observed by the Kepler pho- 
tometer between 1 May 2009 and 15 Jun 2009. 
The upper panel shows the full 44-day time se- 
ries after detrending. The bottom panel shows the 
light curve folded with the orbital period; differ- 
ent plot symbols denote odd- and even-numbered 
transits. The lower curve shows transit data over- 
plotted on the fitted transit model. The upper 
curve covers the expected time of occultation, with 
the fitted (essentially constant) model overplot- 
ted. This model assumes a circular orbit. The full 
folded light curve (not shown here) gives no evi- 
dence for an occultation at any phase. For transit 
and orbital parameters, see Table 1. 
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(|2010l l: iBatalha et aD (|20ld ). Experience with 
both ground- and space-based transit observations 
shows that a fairly large fraction of events that re- 
semble planetary transits are actually caused by 
eclipses involving only stars, or eclipses with prop- 
erties that are significantly confused by diluting 
light from one or more stellar companions, either 
physical or projected. For this reason, we designed 
our follow-up observations to determine whether 
the light curve dips can be ascribed to a transiting 
planet; if not, what eclipsing star or other process 
might be responsible for them; and if so, what the 
properties of the host star and its planet might 
be. 

Ground-based visible-light speckle imaging 
from the WIYN Telescope, NIR adaptive-optics 
imaging from the Mt. Palomar 5m telescope, and 
also with the NIRC2 camera on the Keck tele- 
scope all show Kepler-4 as an isolated star, apart 
from a neighbor that is 3.5 magnitudes fainter at a 
distance of 11.9". We estimate that this star con- 
tributes 2±2% to the flux we measure for Kepler-4, 
which dilutes the transit by a factor of 1.02±0.02. 
We take this dilution and its uncertainty into ac- 
count in computing the the planetary radius Rp 
and its probable errors. Limits o n nearb y com- 
panions are described by Batal ha et al.l (fioiO,); 
the imagery rules out background eclipsing bi- 
naries that might simulate the observed transits, 
up to a magnitude difference of 9.8 in H band, 
for companions between 0.12 and 3.0 arcsec from 
Kepler-4. We have also searched for motion of the 
image centroid that is correlated with the transits, 
finding no such motion with a limit of 4 x 10^^ 
arcsec. 

Two reconnaissance spectra obtained with the 
TRES spectrograph on the 1.5-m Tillinghast Re- 
flector at the Whipple Observatory showed a ve- 
locity variation of less than 150 m/s over five days. 
Accordingly we obtained radial velocity measure- 
ments with the HIRES spect rograph on the Keck 
I telescope (jVogt et al.lll994l '). Figure 2 shows the 
radial velocity data for Kepler-4, folded with the 
transit period and shifted so that transit center 
occurs at zero phase. Assuming zero eccentricity, 
we obtained a radial velocity variation with phase 
that is consistent with the observed light curve, 
with a reflex velocity amplitude K of 9.3^i[ g ms~^ 
and velocity residuals of only 3.6 ms^^. We also 
performed a fit in which we allowed the eccentric- 



ity e to float. This gave e = 0.22 ± 0.08 and K = 
10.0 ms~^; the reduced improved only slightly. 
Finally, a Monte Carlo bootstrap (described be- 
low) that simultaneously fits all photometry and 
RV values, and that is consistent with stellar evo- 
lution models, gives esinw = —0.003 ± 0.058 and 
ecosw = -0.005 ± 0.057. There is thus only 
marginal evidence for a noncircular orbit. In the 
analysis below, we do however account throughout 
for the possible range of e and w in our estimates of 
stellar and planetary parameters and their uncer- 
tainties. Given the integration times of the spec- 
tra and Kepler-4's brightness, T^s, and slow ro- 
tation (wsini = 2.2 ± 1.0 kms~^), the residuals 
are consistent with expectations. We also used 
the HIRES spectra to search for variations in the 
shapes of line bisectors; the component of the bi- 
sector span that is in phase with the orbital period 
has amplitude —1.2 ±4.2 ms~^ (with uncertainty 
estimated from the bisector scatter), offering no 
support for the presence of a blended eclipsing bi- 
nary star. 

We can reject many photometric blend scenar- 
ios, but not all. Kepler-4b's long transit duration 
rules out hierarchical triple systems in which the 
primary of the eclipsing pair contributes less than 
about 25% of the system luminosity. Near-twin 
star systems permitted by this constraint would 
lead to errors in the planet radius and mass by 
factors of a few - these would be serious, but not 
large enough that transits by stars could be con- 
fused with those by planets. Moreover, virtually 
all transiting planet parameter estimates are vul- 
nerable to confusion of this sort. An unresolved 
background eclipsing binary star might accurately 
simulate the transit light curve, but such a sys- 
tem would not likely produce the very small ob- 
served RV variation, nor the small line bisector 
variation. The most plausible remaining blend 
scenario involves an unresolved background tran- 
siting Jupiter-size planet. Such a blend could be 
consistent with all current observations, though 
the required coincidence in apparent position is 
very unlikely a priori. 

Thus, although some kinds of confusing pho- 
tometric blends are possible, there is no evidence 
for them. In what follows, we assume that the ob- 
served light curve results from transits by a rather 
small extrasolar planet across the face of a normal, 
single. Sun-like star. 
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3. Properties of the Planet and Host Star 
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Fig. 2. — Top panel: The phased radial velocity 
curve of Kepler-4, consisting of 19 epochs observed 
using the Keck/HIRES spectrometer, spanning 69 
days. The dashed overplotted curve is a fit as- 
suming a circular orbit, phased to match the tran- 
sit photometry. The solid curve shows the best- 
fit eccentric orbit, with e — 0.22. Middle panel: 
0-C residuals of velocities relative to the circular 
orbit fit. Bottom panel: Bisector span for each 
epoch, measured between the 40% and 90% depth 
points in the cross-correlation against a represen- 
tative spectrum. 



3.1. Method for Estimating Host Star 
Masses and Radii 

The Kepler photometer produces a light curve 
that implies the planetary radius in units of the 
host star radius. Likewise, groundbased followup 
Doppler spectroscopy yields the mass of the planet 
in terms of the host star's mass. Thus, having 
accurate values for the host star's fundamental 
parameters - mass, radius, effective temperature, 
and composition - will be essential to carrying out 
Kepler's mission to characterize planets circling 
distant stars. 

But it is also possible to reverse this logical 
flow, and us e transiting planets as pro b es of their 



host stars. ISeager &: Mallen-OrnelasI (|2003f ) de- 



rived expressions relating the light curve to the 
mean density of a transiting planet's host star. 
They showed that can be expressed in terms 
of the planet's orbital period, the fractional flux 
obstructed by the planet, and two different mea- 
sures of the transit duration. All of these quanti- 
ties are directly measurable from the light curve, 
and their interpretation depends only on Newto- 
nian mechanics and geometry. 

The method we use (henceforth, the p^, method) 
for estimating stellar radii and other parame- 
ters is an im p lemen tation of that described by 
Sozzetti et al.l (120071) and subsequent l y use d by. 



e.g.. iBakos et al.l (l2007f): 



ICharbonneau et al.l ( 2007 ). This technique has 



Winn et al.l (|2007f l and 



become the most trusted way of obtaining infor- 
mation abo ut transitin g planets' host stars, as il- 
lustrated by iTorres et all (|2008l ) . 

The basis for the method is that, although 
stellar models are ordinarily taken to depend 
on 5 parameters {mass, age, metallicity, ini- 
tial helium abundance, and mixing length} = 
{M^,,A^,, [Z],YQ,a}, the last two of these do not 
vary much among stars, and for practical pur- 
poses are often taken as known. In particular, the 
much -used Yonsei-Y ale (henceforth YY) model 
grid (|Yi et al.l 120011 ) uses specified values for the 
mixing length and for the initial helium abun- 
dance. In this approximation, stars are described 
by a 3-dimensional grid of models, parameterized 
by mass, age, and metallicity. Every transiting 
planet discovered by Kepler and followed up with 
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reconnaisance spectroscopy has available a transit 
light curve that constrains p,, and also spectro- 
graphic estimates of T^g and [Z]. Moreover, for 
main-sequence stars, the problem of inferring stel- 
lar structure parameters from the 3 observables 
{p,, Teff, [Z]} is usually well-posed, allowing pre- 
cise conclusions without important degeneracies. 
Thus, the quantities that we may readily observe 
usually suffice to isolate a single set of model pa- 
rameters {mass, age, metallicity}; from these, we 
may compute any other global property (e.g. ra- 
dius, luminosity, or surface gravity). 

We implement the method by searching a 
precomputed grid of models (interpolated from the 
YY models) to find the one that best matches the 
observations in a sense. The optimization prob- 
lem involves interpolating within the given model 
grid and performing a conventional (eg "amoeba" ) 
hill-climbing optimization search. Estimating er- 
rors and verifying that the search has converged to 
the global optimum is done using a Markov Chain 
Monte Carlo (MCMC) process. To allow for pos- 
sible systematic errors in the spectrographic anal- 
ysis, we took uncertainties in [Z] to be the larger 
of 0.06 dex and the quoted value, and twice the 
quoted value for Toff. For p*, we used a param- 
eterization of the (often highly asymmetric) er- 
ror distribution estimated from a jack-knife anal- 
ysis of the photometry. Not all parameter choices 
{M^,^^, [Z]} permitted by a given are consis- 
tent with the photometry, and vice versa. More- 
over, in the light curve analysis, the estimated 
Ph. depends somew hat on the initial guess for AI^, 
(|Koch et al.ll2010h and on the range of orbital ec- 
centricity allowed by the observations. To account 
for the different information supplied by these two 
kinds of model fitting, we iterated the solution by 
putting the MCMC probability distributions back 
into the light curve analysis. The most likely pa- 
rameter values and uncertainty distributions we 
report below are the result of this once-iterated 

fit. 

Brown diooi) has tested the p* method against 



a sample of 169 stars (mostly member s of eclipsing 
binari es taken from the compilation bv lTorres et al 



;hat have accurately known masses and 
radii. This comparison shows that the method's 
systematic errors do not exceed 2-3% in radius, 
and about 6-9% in mass. Random errors ap- 
pear to arise about equally from uncertainties 



in the estimated stellar T^g and metallicity, and 
from uncertainties in the masses and radii in- 
ferred for the eclipsing binary components using 
traditional (but largely model-independent) tech- 
niques. Systematic errors arise mostly because 
rapidly-rotating and hence magnetically active 
stars are seen to have larger radii than slowly- 
rota ting stars wi t h oth er wise similar pr o pertie s 
('eg. [Torres et all (120061 ): lL6pez-Mora"S pOOVl l: 



Morales et al 



"(|2008l )). The sample of eclipsing 
binary stars consists almost entirely of such fast 
rotators; this causes an excess of up to several per- 
cent in the actual radii of stars of sub-solar mass, 
relative to the YY models. On the other hand, the 
host stars of confirmed transiting planets found by 
Kepler are mostly slow rotators (if only because it 
is difficult to measure precise radial velocities of 
rapidly-rotating stars). Errors in these estimates 
should therefore be dominated by measurement 
errors, especially in [Z] and in p, . 

3.2. The Star Kepler-4 

Simultaneous fits to the transits of Kepler-4b 
and to the RV data (allowing nonzero eccentricity) 
yield the ratio of orbital semimajor axis a to stellar 
radius i?*, namely a/i?* = 6.47^q'28, which with 
the orbital period gives p* — 0.50 ± 0.16 gcm"'^. 
This is a low value, roughly 1/3 that of the Sun, 
suggesting that Kepler-4 must have evolved con- 
siderably away from the zero- age main sequence. 
Applying the p* method confirms this conclusion, 
and in fact finds a range of model parameters, all 
of which fit the observations almost equally well, 
but which imply distinct evolutionary states for 
the star. Figure 3 shows evolution tracks for three 
models having masses between 1.13 Mq and 1.28 
Mq. In this mass range, old main-sequence stars 
have small convective cores. Because of efficient 
mixing, these cores suffer hydrogen exhaustion all 
at once, following which the entire star must com- 
press and heat until the central temperature rises 
enough to start shell hydrogen burning. The blue- 
going hook at main-sequence turnoff is a result of 
this compression. In the case of Kepler-4, the ob- 
served p, and Teff (indicated by the error box in 
Figure 3) are closely matched by each of the evolu- 
tion tracks, albeit at different ages and with mod- 
els occupying different evolutionary states. Thus, 
for the lowest-mass star to fit the observations, it 
must have just finished its core-contraction phase, 
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and be on its way to becoming a shell-burning 
subgiant. The highest-mass star to fit the same 
observations must still be on the main sequence 
but nearing hydrogen exhaustion in its core, just 
poised to begin core contraction. 

It is not possible to distinguish among these 
possibilities using existing observations. Although 
the masses of acceptable models differ by as much 
as 13%, the radii differ by only 1/3 as much (in 
order to give the same observed mean density). 
This radius difference would cause a difference in 
log(g) that is too small to measure with current 
methods. Likewise, the luminosity difference of 
about 8% implies a parallax difference of only 4%, 
also too small to discern at this star's likely dis- 
tance of about 550 ±80 pc. Asteroseismology may 
however offer a way out of this quandry. Because 
of the star's relatively large luminosity and ap- 
parent brightness, it may be possible to measure 
the frequencies of its pulsation modes, and hence 
distinguish among the feasible evolution scenarios. 
Kepler-4 is now being observed with Kepler's short 
(60 s) cadence; the mission will report pulsation 
properties estimated from these data if and when 
the pulsation signals rise above the noise. 

We list the inferred properties of Kepler-4 in 
Table 1. For the quantities and R^, we give un- 
certainties that include the effects of observational 
uncertainties, of uncertainty in orbital eccentricity, 
and of uncertainty in the host star's evolutionary 
state. 

3.3. The Planet Kepler-4b 

Given the properties of its host star, the depth 
of transits due to Kepler-4b implies a planetary 
radius of 0.357 ± 0.019i?j = 3.99 ± 0.21 i?©, and 
a mass of 0.077 ± 0.012Afj = 24.5 ± 3.8 M®. As 
with M-i, and i?*, these values and uncertainties 
should be interpreted as describing the centers and 
68%-probability points of the marginal probability 
distributions, accounting for all of the uncertain- 
ties. Kepler-4b is therefore slightly more massive 
than Neptune, and about the same size. Its mean 
density is about 1.9 gcm^"^, greater than Jupiter's 
and considerably larger than Saturn's. Assuming 
a 10% Bond albedo and efficient heat redistribu- 
tion to the planet's nigh t side, its equilibri um tem- 
perature is 1650±200K (|Koch et al.ll2010f ). Radial 
velocity data provide no evidence for other mas- 
sive planets in small orbits. Because of the limited 



time span and precision of the extant RV obser- 
vations, it is however impossible to rule out such 
planetary companions. 

4. Discussion 

Figure 4 shows Kepler-4 in a mass-radius di- 
agram that spans the range so far occupied by 
small transiting exoplanets. Kepler-4b is the 
third known transiting Neptune- like planet, to- 
gether with GJ436b and HAT-P-llb; all three 
have masses and radii equal or larger than Nep- 
tune and Uranus. The most important differences 
between the 3 exoplanets are their host stars and 
the incident flux the planets receive (GJ436b: 
M dwarf, Teq = 650 K; HAT-P-llb: K dwarf, 
880 K; Kepler-4b: G subgiant, 1650 K). On the 
mass-radius diagram Kepler-4b and GJ436b both 
lie below the Z = 0.9 Baraffe et al.(2008) non- 
irradiated model. In these models, Z = 0.9 de- 
notes that the planet has a 10% by mass envelope 
of H and He, and 'non-irradiated' corresponds 
roughly the case of GJ436b and cooler. The fact 
that Kepler-4b and GJ436b have essentially iden- 
tical rad ii (3.88 ± 0.15 i?m fo r GJ436b from NASA 
EPOXI (|Ballard et al.ll2009l ) as compared to 3.99 
± 0.21 i?0 for Kepler-4) at the same mass, de- 
spite Kepler-4b's high Toq, implies a difference in 
bulk composition. We suggest that Kepler-4b has 
a H/He envelope of about 4-6% by mass and a 
correspondingly higher water and rock fraction. 
However, the inherent degeneracy in this part of 
the mass-radius diagram means that the H/He 
envelope might be slightly more massive if the 
rock-to-water ratio in the interior of Kepler 4b is 
unusually high (i.e., significantly higher than in 
Neptune or Uranus). Nevertheless, we can state 
with a measure of confidence that there are no pos- 
sible interior models for Kepler-4b with no H/He 
envelope and neither it nor GJ436b is compact 
enough to be a water-rich super-Earth. 

The slightly eccentric orbits of the two previ- 
ously known transiting Neptune-like planets, in 
similar potentials and at similar ages, suggest that 
their interiors are somewhat less dissipative to 
stellar tides than the giants in our Solar System. 
Therefore, it would be very interesting to resolve 
whether the orbit of Kepler-4b has comparable ec- 
centricity or is circular; the current observations 
are still inconclusive. 
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Table 1 

System Parameters for Kepler-4 



Parameter 



TYansit and orbital parameters 
Orbital period P (d) 
Midtransit time E (HJD) 
Scaled semimajor axis a/Ri, 
Scaled planet radius Rp / Ri, 
Impact parameter 6 = a cos i/R^ 
Orbital inclination i (dcg) 
Orbital semi-amplitude K (ms~^) 
Orbital eccentricity e 
Center-of-mass velocity 7 (ms~^) 

Observed stellar parameters 
Effective temperature Tgff (K) 
Spectroscopic gravity log fir (cgs) 
Mctallicity [Fe/H] 
Projected rotation i; sin i (kms^""") 
Mean radial velocity (kms~^) 

Derived stellar parameters 
Mass M*(M0) 
Radius R^{Rq) 
Surface gravity log 5* (cgs) 
Luminosity (-^o) 
Absolute V magnitude My (mag) 
Age (Gyr) 
Distance (pc) 

Planetary parameters 
Mass Mp (Mj) 
Radius Rp {Rj, equatorial) 
Density pp (gcm~^) 
Surface gravity log^p (cgs) 
Orbital semimajor axisa (AU) 
Equilibrium temperature Teq (K) 



Value Notes 



3.21346 ±0.00022 A 

2454956.6127 ± 0.0015 A 

"■^'-0.28 

02470+""""^^ A 

U.UZ4t(U_Q QQQ3Q A 

om2tl-l%\ A 

9.3;j;^ A,B 

(adopted) A,B 

-1.27 ±1.1 A,B 

5857 ±120 C 

4.25 ±0.10 C 

+0.17 ±0.06 C 

2.2 ±1.0 C 

-61.0 ±0.10 B 

1.487ig;SJ^ C,D 

4.17 ±0.04 C,D 

2.26j:°;^| C,D 

4.00 ±0.28 D 

4.5 ± 1.5 C,D 

550 ±80 D 

0.077 ±0.012 A.B,C,D 

0.357 ±0.019 A.B,C,D 

1.91^°;^? A,B,C,D 

3.1612'Jo A,B,C,D 

0.0456 ± 0.0009 E 

1650 ± 200 F 



Note. — 

A: Based primarily on the photometry. 
B: Based on the radial velocities. 

C: Based on spectrum analysis (FIES/MOOG or HIRES/SME). 

D: Based on the Yale-Yonsci evolution tracks. 

E: Based on Newton's version of Kepler's Third Law. 

F: Assumes Bond albedo — 0.1 and complete redistribution. 
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6000 5900 
Teff (K) 



Fig. 3. — Yonsei-Yale evolutionary tracks for stars 
with [Z] = 0.17 and masses between 1.13 and 1.28 
Mq, plotted in the — Teff plane. The box in- 
dicates our adopted observational constraints on 
and Teff. For this range of masses, all tracks 
lie partly inside the error box, showing that mul- 
tiple solutions are possible that satisfy these con- 
straints. 



Baraffe Z=0.5, NI, TGjr _ 



■ HD 14902^ 



Baraffe Z=0.9. 1, 3,5,7 Gyr: 




_ Baraffe Z=0.9, NI. 3,5,7 Gyr : 

HAT-P-llb : 

[^Ke pler-4b 

' Uranus □ 'i-BT ~ 
^^P'^""^ [ ValenciaZ=l,507<.water ■ 

GJ436b 



2 -- CoRoT-7b 



Valencia Z=l, rock 



20 30 40 50 60 70 80 90 100 

Plaiietai-y Mass (Earlli Masses) 



Fig. 4. — Kepler-4b is shown in a mass- 
radius plot, along with other hot Neptunes 
the hot super-Earth s CoR oT-7b (jLeger et al 



20091: iQueloz et al 



20091: 
(|Cha r 

Valeni 



Charbonneau et al 



.20091) and GJ1214b 
20091). and model curves 
fc BarmanI tOO^ and by 



^ arafFc. Cha rbrier. .. 

Valencia. Sasselov. fc O'Connelll tOO'K showing 
increasing heavy-element fractions toward the 
bottom of the Figure. Multiple dashed and solid 
curves show model results for differing ages; 
"I" indicates irradiated and "NI" non-irradiated 
models. Kepler-4b appears to be denser than 
HAT-P-llb but similar to GJ 436b, denser 
than the non-irradiated Z = 0.9 model by 
Baraffe. Charbrier. fc BarmanI (|2008r ). but much 
less dense than is expected for a water or rocky 
planet. 
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